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Non-positive Markov approximations are sometimes used to describe the dynamics 
. of qubits in weak interaction with suitable environments; the appearance of negative 

probabihties is avoided by assuming that the transient regime ehminates from the 
possible initial conditions those qubit states which would otherwise be mapped out 
of the Bloch sphere by the subsequent Markovian time-evolution. By means of a 



O . 

\^ , simple model, we discuss some physical inconsistencies of this approach in relation 

o , 

, to entanglement; in particular, we show that slipped non-positive reduced dynamics 



might create entanglement through a purely local action. 



Cr. I- INTRODUCTION 



^ ' Semigroups of dynamical maps are used to describe the time-evolution of open quan- 

tum systems S in weak interaction with suitable external environments, typically an 
infinite heath bath in equilibrium at a given temperature acting as a source of dissipation 



nenomeno 



mm 



ogical applications in 



and noise. They have been successfully used in many p 
quantum chemistry, quantum optics, statistical physics 

In the following, we shall consider qubits described by density matrices p, correspond- 
ing to three dimensional real vectors of length < 1 in the Bloch sphere, that evolve in 
time under the action of semigroups of linear maps 7t, t > 0. Formally, these arise 
from the exponentiation of a generator L, 7^ = exp(tL), and satisfy the forward in time 
composition law 74 o 7^ = •jt+s, t,s >0. 



A preliminary natural request on the maps 74 is that they send any initial state p into 
another state 7t[p] at all t > 0, namely that they map the Bloch sphere into itself. Only 
in such a way, the spectrum of an evolving p remains positive and its eigenvalues can be 
interpreted as probabilities. 

Such a property of the maps 7^ is called positivity. In line of principle, it is not sufficient 



to guarantee full physical consistency o 

nnn 

complete positivity, need be imposed p, 121 



the maps 74: a more restrictive property, namely 
In such a way, not only the positivity of jt is 



guaranteed, but also that of the amplified map 7*® id; this describes the time-evolution of 
a qubit S statistically coupled to an ancillary qubit which remains inert under the action 
of the identity operation id. Moreover, complete positivity of 7^ fully characterizes the 
form of the generator L P, 0| . In particular, in the case of a qubit system, there appears 
a characteristic order relation (2Ti > T2) between the decay times of the diagonal (Ti) 
and off-diagonal (T2) entries of its time-evolving density matrix j^. 

The semigroup 74 describes the reduced dynamics of the immersed qubit S when the 
environment degrees of freedom have been eliminated and the memory effects due to a 
short transient regime have been got rid of by means of suitable Markovian approxima- 
tions. If not performed with due care, these latter lead to reduced dynamics neither 
completely positive, nor even positive. 

Absence of positivity of jt means that there are initial density matrices p that may, in 
the course of time, develop negative eigenvalues and thus lose their meaning as physical 
states. If one wants to stick to non-positive reduced dynamics, a possible way out of 
physical inconsistencies amounts to assuming that not all density matrices are allowed 
as initial conditions for 7f, but only those which do not develop negative eigenvalues. 
The mechanism which eliminates the unwanted initial states is ascribed to the transient 
regime which rules the time behavior of the subsystem S before one can legitimately use 
the semigroup 7^. Namely, prior to 7^, the environment action on the subsystem is via 
a map S that projects the whole state space of S, S{S), into a subset of "good" states 
§(iS(S')), so that 7f o § acts as a positive map on S{S), even if 7^ does not. 

The map S is known in the literature as slippage of initial conditions jll . 
its introduction being motivated by the the difficulty to accept that a physical effect like 
the decay-times hierarchy of a qubit be induced by its possible entanglement with an 
uncontrollable external ancillary qubit, a seemingly academic, abstract scenario with 
philosophical overtones. 
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On the other hand, quantum information and communication theories have amply 
demonstrated the role of entanglement as a concrete physical resource and developed 
techniques both theoretical and experimental for its manipulation Q|. It is thus not only 
of academical interest to study the slippage- approach in relation to entanglement; more 
precisely, we shall be interested in setting the ground to answer the following question. 
Suppose the inconsistencies of a non-positive 7^ have been cured by considering a slipped 
7t o §. Which kind of effects has (74 o §) (g) id on the entangled states of S + S7 

Aim of this paper is to extend a previous result 12| in order to indicate which kind of 
pathologies may arise from amplifications of slipped non-positive semigroups acting on 
bipartite systems. Given an environment, its associated slippage operator § is far from 
being technically accessible as a mathematical object and such are its effective properties. 
We shall thus model § by the simplest map with the prescribed slippage properties and 
then discuss its action on the isotropic states [l^ of a two qubit systems, one qubit 
immersed in a dephasing environment and the other one external to it. We shall show 
that, by preceding it by a suitable slippage operator §, a non-positive reduced dynamical 
map 7t can indeed be turned into a completely positive map 7^ o S, but also that, unless 
strongly slipped, the maps 74 can, acting locally, create entanglement, a fundamentally 
non-local property. Furthermore, we shall present instances of cases where the necessary 
strong slippage results in the elimination of all isotropic entangled states. We take this as 
an indication of a general behavior: in the weak-coupling regime, slipped non-completely 
positive dynamics appear to be incompatible with entanglement. 



II. QUANTUM DYNAMICAL SEMIGROUPS 

In this section, we shall briefly review the standard approach to open quantum dy- 
namics, namely the so-called weak-coupling limit. 

The typical physical context is represented by a finite-dimensional (ra-level) subsystem 
S in weak interaction with an environment E, this latter corresponding to an infinite 
dimensional reservoir in equilibrium at some fixed temperature T. The space of states 
S{S) of S consists of n x n density matrices ps, while the state of the environment pe 
is taken to be an equilibrium state pe- The dynamics of the (closed) system S + E is 
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typically described by a Hamiltonian of the form 

Hs+e = Hs + He+ gHi , = Y,X''s ® X% , (1) 

a 

where Hs is the subsystem Hamiltonian, Hi an interaction term, of strength g, linear in 
the system and (centered) environment operators Xg^ {Tie{.PeX%) = 0). 

When the interaction between S and E is weak {g << 1), it makes sense to derive 
a so-called reduced dynamics, namely a description of the time-evolution of S involving 
the system 5* alone. The derivation amounts to the elimination of the degrees of freedom 
of E while retaining the effects of their presence on S; the resulting time-evolution is 
irreversible, in general non-linear, and dominated by memory effects. By assuming the 
initial state of S" + to be of the uncorrected form ps® Pei non-linearities are automat- 
ically ehminated and the standard projection technique leads to the following equation 
of motion for ps e S{S) Q: 

\g' [\sK{s)[ps{t-s)], (2) 



dtpsit) 



Hs,Ps{t) 







where K(s) is a highly complicated kernel that, acting on the states of S, takes into 
account both the degrees of freedom of E and the time-evolution of all of them prior to 
time t. 

The ensuing dynamics preserves the trace of ps and is completely positive in the sense 
that it sends any initial ps into a ps{t) which results from the action of a linear map Gf- 

p ^ psit) =: = J2^jit)psV;it) , (3) 

j 

where the Vj{t) are n x n matrices such that J2j ^^(^)^ (^) = 1- 
Remarks 2.1 

1. A trace-preserving linear map G on S{S) is called positive if G[p] has positive 
spectrum for any p G S{S). If G is completely positive, then it is not only positive, 
but, upon coupling S with another n-level system system, the amplified map G(S>id 
preserves the positivity of all possible states of the compound system S + S, where 
id means that the ancillary system S is not affected. It turns out that complete 
positivity is equivalent to G being of the Kraus-Stinespring form [3, ^ 

G[p] = J2G,pG], J2^]G, = 1. (4) 







2. The physical meaning of complete positivity is intimately related to the phe- 



nomenon of quantum entanglement p; therefore, its physical role can be appreci- 
ated only considering correlated bipartite quantum systems of which only one party 
undergoes a state change. In fact, if coupling of S to ancillas could be excluded, 
the physical transformations of the states of S might adequately be described by 
positive linear maps G. Also, if S is coupled to another S, amplified positive maps 
G ® id would preserve the positivity of any separable state 

PZs■■=Y.^^JPs®P's, P^^^iS). (5) 

However, since there are states of S+S, the entangled ones, which cannot be written 
as in (0), when G (8> id acts on anyone of them, the positivity of the eigenvalues 
of G ® id[p5"(^5] can be preserved if and only if G is not only positive, but also 
completely positive. 

3. If a linear map G on the states of S is not completely positive, this does not mean 
that G ® id sends any entangled state of 5* + 5* out of the class of states: some 
entangled states may none the less have the positivity of their spectrum preserved 
at all times. However, for any non-completely positive G, the totally symmetric 
projection 

1 " 

^■=-Yl ' (6) 

will always have some negative eigenvalues appearing in the spectrum of 7^ (S> id[P]. 



Indeed, a theorem of Jamolkiowski 



18| establishes that G ® id[P] is a positive 



matrix if and only if G is completely positive. 



Beside being practically impossible to handle mathematically, the dynamical maps Gf 
do not satisfy a forward-in-time composition law: GtoG^ 7^ Gt+s- Manageable Markovian 
approximations are obtained by weak-coupling limit techniques 2!; these exploit the 
weakness of the interaction between S and E. Indeed, (j21) implies that the influence of 
the environment becomes visible on a time scale such that g^t = r; by setting t = r/g"^ 
in (j21), the typical prescription is to let (7 — and approximate by and p{t — s) = 
pir/g"^ — s) by p{r j g^^ = p(t). One thus obtains the equations 



dtPsit) 



Hs,Psit) 



+ g^]D)[psit)] , 



(7) 



D 



where D := dsK(s) and Hs is a Lamb-shifted subsystem Hamihonian corrected by 
terms of order g'^. 

The hnear map D takes into account the environment induced dissipative effects, 
essentially damping and noise, affecting the subsystem S when it weakly interacts with 
its environment E. By introducing a Hilbert-Schmidt orthonormal set of n x n matrices 
Fj, i = 1, 2, . . . , — 1, = l/y/n, Tt^F-Fj) = D can always be put in the form 

B[ps] = C,,[f,psF} - ^{f]F,.p 



where C = [Cij] is a self-adjoint (n^ — 1) x (ra^ — 1) coefficient matrix, the so-called 
Kossakowski matrix [9j. 

By setting the right hand side of (jHl) equal to L[p5(t)], the resulting Markovian reduced 
dynamics consists of a semigroup of linear maps 74 = exp(tL); these maps are trace- 
preserving and their positivity or complete positivity depends on the matrix C: 

1. if the Kossakowski matrix C > 0, then the maps 74 are completely positive and 



vice versa 



221; 



2. if C is not positive, nogeneral necessary and sufficient conditions are known which 
yield positive 7t (see p for certain cases where this possibility indeed exists). 



III. AN OPEN QUBIT SYSTEM 

We now follow and develop an argument developed in and consider, as a con- 
crete model, a single qubit, represented by a spin 1/2 particle coupled to a constant 
magnetic field vertically directed and to a weak stochastic classical magnetic field 
G = {Gi(t),G2{t),G3(t)). The time-evolution of the qubit density matrices is deter- 
mined by the time-dependent Liouville-von Neumann equation 



dtPsit) 



^^0-3 , Ps{t) 



— I 



Gi{t) (Ji , ps{t) 



i=l 



(9) 



where the Uj are the Pauli matrices and the scalar quantities Gi{t) are chosen to be 
stationary, stochastic variables with vanishing first moments and diagonal covariance 
matrix: 

{Gi{t)G^{s)) = G.^-^^^'-'^5i, , (10) 



with Gi > G2 > 0, 6*3 > and Ai = A2 = A > 0. Such a classical stochastic field is a 
convenient modeling of a heat bath whose temperature is sufficiently high, whereby the 
global Hamiltonian for the system plus environment may be replaced by Q. 

Of course, the solution of the above stochastic equation is a stochastic density matrix; 
useful physical information may only come by considering the average of ps{t) with 
respect to the noise. Techniques based on the weak-coupling limit can be applied to (jH)) 
(l9| ; in general, these lead to a reduced dynamics which suffers from lack of positivity of 
the ensuing Markovian reduced dynamics. 



With reference to the evolution equation (|7j) with the dissipative term as in (jH}, choos- 
ing Fi = i = 1,2, 3, and Fa = —;=, one gets 

H = u;as, uj:=u{1+^-^^1±^\ (11) 



C 



^ ai -b 0^ 

—b a2 
y ay 



A2 + 4^52 



2G12X 2G3 „~ G2 — Gi 



In order to better expose the consequences of sticking to non-positive dissipative semi- 
groups, we simplify the analysis by choosing u; ^ A3 ^ A and XGi^2 ^ i^\G2 — Gi\. 
This allows us to neglect 0:1,2 with respect to both b and a; we shall thus deal with the 
dissipative evolution equation 

dtPs{t) = -i[uja3, ps{t)] + a(^a3 ps{t) - ps{t)^ -b (^(Ti ps{t) a2 + cr2 ps{t) . (13) 

In order to solve the above equation and study the ensuing time-evolution, it proves 
convenient to recast the qubit density matrices as 

= lil + r-a), (14) 

where a = (ai, 0-2, 0-3) and r = (ri, r2, r^), with 

n = 2^{pu) , r2 = -2^ipu) , r3 = 2pn - 1 • (15) 
The matrix p is positive if and only if its determinant 

Det(p) = > , (16) 
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whence each state is identified by a vector r G of norm ||r|| < 1, that is by a point in 
the so called Bloch sphere. 

If the qubit is weakly interacting with an environment, its dissipative dynamics should 
be described by a semigroup of trace-preserving maps jt sending any state p at t = 
into a state '~ft[p] at time t > 0. According to ()16|) . this is equivalent to \rt\ < 1, where 

= {f^i{t),i^2it),r3{t)) G is the Bloch vector identifying 

1 ^ 1 ^ 

^t[p] = ^{l + J2'^^ltW^]) = -(l + $^r,(t)a,) . (17) 



i=l i=l 

In particular, ()13|) translates into a linear differential equation for the Bloch vector 



-2C 



(18) 



where the 3x3 matrix C = Ti + T> consists of the sum of an antisymmetric component 
Ti corresponding to the commutator with the Hamiltonian in ((7j), while T) is symmetric 
and corresponds to the dissipative term D. Explicitly, 



^ a h + uj ^ 



b — to a 




n 



uj 
-cj 




a b 
b a 
yO Oy 



(19) 



The semigroup of maps jt = e*'^ on the state-space 15(5") generated by (fT!^ corresponds 
to a semigroup of 3 x 3 matrices Qt = e~^*^ acting on the Bloch sphere. 
From (UHl), (Uni) and (HSl) it follows that 



Ti = —2a ri — 2{b + u) r2 
r2 = —2(6 — Lj)ri — 2a r2 ■ 
r3 = 



(20) 



Setting Q'^ = uj'^ — 6^ one gets: 



nit) -- 

r2{t) = e" 



-2at 



2at 



ri cos(2^]^) - r2^ sin(2fit) 



sin(2fit) + cos(2fit) 



(21) 



and 



Apparently, the action of the semigroup is that of a dephasing channel with no in- 
fluence on the diagonal elements of any initial p: it remains to be checked whether the 
matrices Qt correspond to non-positive, positive or completely positive maps 7^. 

The Kossakowski matrix in (fT^ is positive, hence 7t completely positive, if and only 
if 6 = 0. 

Let instead 6 7^ and consider the Bloch vectors y± = l/v^(0, 1, ±1, 0) corresponding 
to the pure states (Det(p±) = 0) p± = |(1 + o"i ± (T2)- At t > they evolve into matrices 
of trace 1 with Bloch vectors such that 



|r±||2 _ -Aat 
II ~ ^ 



cos2(2fit) + ^ sin2(2fit) t tt sin(4^]^) 



U{a ± h) , (22) 



for t 0. It thus follows that 7^ acts in a physically consistent way on the states of the 
qubit only if — 6^ > 0. Actually, this condition is also sufficient for 7^ to preserve the 
positivity of states for all t > 0. Indeed, — 6^ > corresponds to the positivity of the 
matrix V in ()19|). Now, a state p exits the Bloch sphere at time t if and only if the Bloch 
vector of 7t[p] is such that = 1 and 

= -2{vt\V\vt) = -2a{rl{t) + rl{t)) - 2hr^{t)r2{t) > . (23) 

Thus, when V > 0, the state moves from the surface towards the interior of the Bloch 
sphere if {rt\V\rt) > 0, while if {rt\V\rt) = 0, then it remains of norm 1. 

Proposition 3.1 The semigroup Qt = e~^*^ generated by (fTHj) with C as in p9|) corre- 
sponds to a semigroup of linear maps 7* on the state-space of the qubit 5* which are 

1. completely positive if and only if 6 = 0; 

2. only positive if and only if & 7^ and — 6^ > 0; 

3. not even positive if < a < 6. 



IV. SLIPPAGE OF INITIAL CONDITIONS 



From the considerations in Remark 2.1.2, in order to avoid physical inconsistencies 
that may arise from non-completely positive reduced dynamics in presence of entangled 
states, it would be rather natural to consider only reduced dynamics enjoying complete 
positivity and to exclude those that do not. Actually, the problem is more acute since 



iU 



those Markovian approximations 7t which fail to be completely positive are sometimes 
not even positive. This means that their pathological behaviour, of the kind discussed in 
the previous section, already appears with a single qubit. 

At the level of a single system, if one wants to stick to non-positive reduced dynamics 
as the best Markovian approximations to the actual dissipative dynamics, a possible way 



out of the p 



ivsica. 



mm 



inconsistencies is known in the literaure as slippage of initial conditions 



The argument goes roughly as follows. Not all density matrices ps evolve into •ytlps] 
with negative eigenvalues in their spectra; thus, let S*{S) C S{S) denote the subset of 
states which remain positive under 7^ for all t > and let § : S{S) 1-^ S*{S) be a map 
selecting this subset of "good states" for 7^. Then, even though the maps •jt suffer from 
physical inconsistencies, the maps 74 o S do not. 

The slippage of initial conditions is supposed to be operated by the transient regime 
prior to the Markov one. In other words, according to 

ptfrans 

§H ^ / / dsK{s)[ps{ttrans - S)] , (24) 

Jo 

where ttrans is the span of time during which memory effects cannot be disregarded. 
Notice that the weak-coupling limit g'^t = t, g —>■ 0, does exactly push tfrans ^ 0. 
Roughly speaking, the action of the transient regime is supposed to be such that the 
Markovian dynamics following it does not act on any possible initial density matrix, 
but only on those selected by S which are free from the inconsistencies arising from 
non-positivity. 

Of course the plausibility of such an approach has to be tested against actual physical 
open quantum contexts; however, the highly complex nature of the kernel K(t) describing 
the transient, memory-full dynamics, scarcely allows any hope to go beyond extremely 
rude approximations in describing mathematically what its effects really are. In the 
following, relative to the reduced dynamics in ()2H1 . we shall model a slippage operator 
in a simple way that however embodies some of its salient features. 

We first observe that according to the slippage philosophy, there would be no need 
to invoke a selection map S when > 6^ for in this case the maps 7^ are positive and 
cannot drive any initial state outside the Bloch sphere. In other words, all initial 2x2 
density matrices are good initial states for the reduced dynamics and one need not resort 
to a slippage mechanism to select a suitable subclass of them. 



ii 

We shall then set o? < 6^; from the arguments of the previous section, the most 
dangerous states are pure, so that non-pure states within the Bloch sphere are expected 
to offer less problems. How far into the Bloch sphere should one go in order to be sure 
to get rid of all those states which 7^ would sooner or later map outside it? It is evident 
from (PT|) that Bloch vectors may exit the Bloch sphere but their lengths ||rt|| cannot 
diverge. The maximum value R achievable by ||rt|| can be explicitly obtained. 

First, observe that, within and on the Bloch sphere. 



t\\ , 



(25) 



where Qf is the transposed of Qt and the operator norm ||^^^t|| equals the largest eigen- 
value of the positive, real symmetric matrix 



+ e 



-4at 



^0 0^ 


+ e-'°'^sin(2m) 



fl o\ 

1 
yO Oy 

^{b-uj) sin(2tn) cos(2tf]) 0^ 



\ 



-n cos(2tr]) {b + uj) sin(2tfi) 




(26) 



Of the three eigenvalues of the above matrix, one is 1 and the largest of the remaining 
two is 



R^(t) := e-^-^* |sin(2tn)| + ^l + ^sm\2m) 
If a'^ < b"^, it can be explicitly computed that R{t) achieves its maximum 



(27) 



R := e 



-2at' 



' UJ 



UJ 



- - 



1 n b^-a^ 

at ^-^arcsm -J^^^ 



:28i 



Remarks 4.1 

1. i? is surely > 1; indeed, for < 6^, 7^ is not positive, whence the norm of Gj^Gt 
cannot be < 1. 



2. The maximum (j^Hj) is achieved at t' starting from the pure state p* whose Bloch 
vector r* is the eigenvector of Q'fQt relative to the eigenvalue R^. 



Set < fi < R ^ and consider the following linear map on S{S): 



P^§M[p] = ^(l + /^r-a) (29) 



Then, because of linearity, 



7t o §^[p] = ^(l + /iFi ■ a) , and ||/irt|| < /xi? < 1 , (30) 

whence no state can be mapped out of the Bloch sphere by 7t o §^ at any t > 0. 
The map is completely positive, its Kraus-Stinespring form being 

3 

Sm[p] = ^P + ^E^^/'^- (31) 

i=l 



Remark 4.2 If the interaction of the subsystem with the environment in equilibrium 
is weak, then the hypothesis of an initial factorized state ps ® Pe is not an unphysical 
restriction, as well as the hypotheisis that the transient act as a linear completely positive 
operator. 

The action of §^ is rather drastic in that it rigidly maps the unit Bloch sphere onto 
a sphere of smaller radius; one could instead envisage more sophisticated slippage mech- 
anisms. However, the intent of this paper is to show the kind of problems afflicting the 
simplest possible of them. They should indeed not depend on the explicit model of tran- 
sient, being it related to the non-positivity of the Markovian approximation and, as we 
shall presently see, to the existence of entangled states. 



V. SLIPPAGE VS ENTANGLEMENT 



Thanks to the rapid development of quantum information, the notion of entanglement 
left the purely epistemological arena in which it was confined for thirty or so years and 
made its appearance on the foreground of physics as a most precious resource to perform 
otherwise impossible tasks as dense quantum coding and quantum teleportation jioj l. 

Physically speaking, the model we are going to consider is that of one qubit immersed 
in a stochastic environment whose effects are described by the maps 7^ of the previous 
section, which is correlated to another qubit which is a dynamically independent and 
inert ancilla. This is the physical context when, for instance, an entangled Bell state is 
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constructed in a laboratory and, while one of the two qubit is kept there, the other is sent 
to a distant party through a noisy channel. The time-evolution of the compound system 
is thus described by the semigroup of maps 7t ® id, where the presence of the bath is 
felt locally by one of the two qubits only and the other does not evolve at all (this is the 
meaning of id). 

In the following, we shall study how the entanglement content of special qubit states 
changes under the local time-evolution as above. Concretely, we shall consider the time- 
evolution of the concurrence [2^ , which, for any state p of a two qubit system, is defined 



where the Aj's are the positive square roots of the (positive) eigenvalues of the matrix 



where p* is the matrix obtained by taking the complex conjugate of the entries of p in a 
chosen representation. 

By means of the techniques sketched above, we shall now study whether curing the 
non-positivity of reduced dynamics by slipping the initial conditions can nevertheless 
conflict with the existence of entanglement between the immersed open qubit and its 
inert ancilla. 

Remark 5.1 That problems may easily arise can be argued by observing that in the 
case of a? > 1)^ > 0, the maps 74 are positive and, a priori, no slippage is needed. 
Nevertheless, according to Remark 2.1.1, the totally symmetric projector P is definitely 
mapped out of the state of space of the compound system S + S, for 74 is completely 
positive only if 6 = 0. Therefore, if 6 7^ 0, one needs a slippage operator act on the two 
qubit system; as the latter can only be due to the bath, it must also affect the single 
qubit immersed in it, independently of whether it is coupled or not to an ancillary qubit. 

In order to set the ground for discussing the slippage mechanism in relation to entan- 
glement, we rewrite the totally symmetric projector as follows 



by 



C(p) := max{0, Ai - A2 - A3 - A4} , 



(32) 



p p , p := ((T2 ® 0-2) p* (0-2 ® (^2) , 



(33) 




1 ® 1 + (Ji (g> (Ji - (72 ® (72 + (T3 ® 0-3 



) 



(34) 
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SO that (jHlll maps it into 

Pf, := §^ O id[P] = ^ ® 1 + yu(^ai ® ai - (72 ® (^2 + O (Tg)^ . (35) 

Then we let •jt ® id act on and consider := ® id[P/^]; the result is easily 

computed from 



-2at 



ai{t) = e 

o-2(t)=e-2-* 

0-3 (i) = 0-3 , 
which follow using ^T^i and (fTTjl . Explicitly, 

/l+yU 2flBt\ 

1 - ^ 2^ a 


\2/iP; i + /xy 



(Ti cos(2f]t) + cT2%^ sm{2nt) 



ai^sm{2Qt) + aa cos(2f]t) 



(36) 



t ■- 



-2at 



cos(2fit) - if, sin(21]t) 



:=i Ag-^"^* sin(2fit) 



The eigenvalues of this 4x4 matrix are, in decreasing order. 



e5'(t) := \ 



em ■■= - 



1 + /i I 1 + 26-^"' \l 1 + ^ sin2(2(]t) 



1 + ^ I 1 - 2e-2'^* ^/ 1 + ^ sin2 (2^]^) 



em 
em 



/i 1 - 2e 



-2at 



n 



sin(2fit) ) 



1 -/i 1 + 2e 



-2at 



- sm{2nt) j 



(37) 

(38) 

(39) 

(40) 
(41) 



As already stressed, physical consistency asks that the slippage operator S^j cure not 
only the non-positivity of 74, but also ensure the positivity of the maps (74 o §^) (g) id 
acting on the the entangled states of the compound system S + S. In particular, P^(t) 
must correspond to a density matrix at all times t > 0. 



Proposition 5.2 Set 



R^{t) := 1 + 2e-^''*- sin(2fit) 



the eigenvalues of P^(t) are positive if and only if 

1 



< < 



Ra 



where Ra 



1 + 2e 



-2aU 



(42) 



(43) 



is the maximum of R4{t) which is achieved a.t t^: = ^ arcsin 



n 



Proof: Since the eigenvalues 63 4 (t) interchange when sin(2i7t) changes sign, one can 



always consider < t < Then, the positivity of 64 (t) requires fiRiit) < 1 for all 
t > 0. Vice versa, suppose /i < -Rj^(t) for alH > and denote by 



the term multiplying fi in the largest eigenvalue e^(t). Since < yU < 1 and 6 > 0, one 
estimates: fi Ri{t) < 2 + Ri{t) < 3, whence 1 > e'^it) > e^(t) > e^{t) > e^(t) > 0. ■ 

In figures 1-3, the red line shows the time behaviour of the largest eigenvalue e^(t) and 
the blue line the lowest eigenvalue e'^{t), for given values of a and b, with /i = 1, /i > R^^ 
and fj. < R^^. Both functions are computed with rescaled parameters a ^ a/to, b ^ b/to 
and if: ut; in this way fi^ = 1 — 6^ with < a, 6 < 1. The green line corresponding to 
height 1 is showed for the sake of comparison with e^(t). 




(44) 
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Figure 1: a = 0.1, b 



0.9, R^ ^ = 0.25, fi = l. Red line: e^(t). Blue line: e^(t) 
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Figure 2: a = 0.1, 6 



0.9, R^^ = 0.25, ^ = 0.4. Red line: e^(t). Blue line: e^(t) 
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Figure 3: a = 0.1, b = 0.9, R^^ = 0.25, /i = 0.2. Red line: e^(t). Blue line: e^(t) 



Remarks 5.2 

1. The states that are obtained from the totally sjTiimetric projector P by the 
action of the slippage operator §^ are the class of isotropic states, P^ = + // P 
in standard form. They are entangled if and only if 1 > > | (l8l |. 

2. According to Remark 2.1.2, condition PHj) makes -P^(t) = (7toS^) (g)id[P] a positive 
matrix for all t > 0, whence all slipped maps 7t o result completely positive, 
despite 7t being possibly not even positive. 

3. If < 6^, we have already seen that the channel can cure non-positivity by 
multiplying the Pauli matrices hj fi < R^^ where R is given in ()28|). From (j^Tj) 
and ()42|1 it follows that R{t) < RA{t) so that R < -R4; namely, the maximal radius 
of the sphere of slipped initial conditions suggested by the single qubit non-positive 
dynamics is too large for keeping physical consistency against entanglement. The 
radius is to be decreased to at least R^^. 

In the following we shall see that even the value does not keep the non-positive 
dynamics free from pathologies; however, the ones we are going to expose are much more 
intriguing and subtler than the appearance of negative eigenvalues, being instead related 
to the creation of entanglement by means of the local action of maps of the form 7^ o 

We shall now consider /i < R^^ which ensures that the -P^(t) are physical states of the 
compound system 5-1- S* at all times t > 0. In order to study their entanglement content, 



we construct their concurrence C^(t) (see (pT^ - (p!?jl ): one checks that -P^(t) = Pfj.it), 
whence the square roots of the eigenvalues of P^(t)P^(t) are the values (PH|) - (PT|) and 



C^{t) = max{0,c^(t)} , c^{t) := ^e"^'^* ^/ l + _ sin^(2fit) ^ . (45) 



At t = 0, := — ^ — > if and only if /i > 1/3, that is if and only if P^ is entangled. 



Proposition 5.3 The states -P^(t) are entangled if and only if 



Rl{t) ^ i?4 ^ ^ 



Proof: The upper bound on fi guarantees that the Pf^it) are well defined density ma- 
trices of the bipartite system 5* + S* at all t > 0, while the lower bound coming from ()46p . 
sets the range of t > for which C^(t) = c^(t) > 0. ■ 

Remark 5.3 Notice that by partially transposing P^(t) in ()37|). one exchanges Bt and 
Ct so that the eigenvalues of T C>?> id[P^(t)], where T denotes transposition, are the same 
of those Pix{t), but with yU — > — /x. Therefore, one checks that corresponds to a non- 
positive partial transposed of P^iit), thus to an entangled P^iit), while C^(t) quantifies the 
amount of entanglement it possesses. 

We are now interested in the change of entanglement with time; from Proposition 5.1, 
we know C^(t) cannot increase if •jt is completely positive. Indeed, 7to§^, the composition 
of two completely positive maps, would also be completely positive and hence a physically 
consistent local operation. 

Let us consider the time- derivative 



2/xe" 



-2at 



c,{t)= ^ G{t), G{t): 

1 + ^sm\2nt) 



bW^^ + a? 



cos(2f2t + if) sin(2i7t) — a 



^]2 

(47) 



where cosv? = /J^ 2 • checks that the function within square brackets achieves its 



maximum 

62 



G := maxG'(t) = — - Vl^M^ - a -a , (45 



ctti t — 1 2 . 



Further, G > if and only if o? < 



This is only possible if < 6^ for cj^ > 6^ 



hence only if 7t is non-positive; in such a case c^(t) increases in a neighborhood of t and 
does it independently of the value of /i. Therefore, if in the same neighborhood > 0, 
then the entanglement of Pfj,{t) increases by the local action of the non positive 74, an 
apparent physical inconsistency, for, as already noticed, physically, entanglement can be 
created by non-local operations only. 

In order to check this possibility, we compare the bound fiR^ < 1, which is necessary 
to physical consistency of the -P^(t) as states at all times, with the lower bound which 
ensures positive concurrence. Since we are interested in a neighborhood of t, we shall 
consider < t + 1 as temporal parameter; then ()4fij) gives 

b 



Ri{t + t) > f{t) := e 



-2at 



52 

1 + — sin2(2(](t + t)) - e 



-2at 



> . (49) 



In Figures 4-7, with the parameters a,b,u rescaled as in figures 1-3, the red lines 
show where f(t) > 0, the blue ones where g(t) := G(t + 1) > and thus the derivative 
in (jlTj) is positive. It is then apparent that there are values of a,b with a < b'^/{2uj) 
(a < 6^/2 in the rescaled parameters) such that there exist time-intervals t G [^1,^2] 
where fl49|l holds together with c^(t + t) > 0. For such choices of a and 6, in order to 
avoid the unphysical creation of entanglement by means of local operations, one must 
exclude /i G + t)~^, R^^] with t G [^1,^2] and must thus enforce the stronger bound 

1 



< /i < 



(50) 



maxte[tj,t2] Ri{t + t) 

Furthermore, in the same figures 4-7, the green lines plot the function r{t) := i?i(t+t) — 3; 
then, figures 5 and 7 show that, in the interval [^1,^2] where f{t) > and g{t) > 0, by 
decreasing a, it also holds r(t) > 0, whence that Ri(i + t)^^ < 1/3 and, consequently, 
also the upper bound in is samller then 1/3. According to Remark 5.2.1, it turns out 
that, in such cases, in order to ensure physical consistency, must destroy all entangled 
isotropic states. 
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Figure 4: a = 0.3, b = 0.8. Red line: f{t). Blue line: g{t). Green line: r(t) 
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Figure 5: a = 0.1, b = 0.8. Red line: f{t). Blue line: Green line: r{t) 
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Figure 6: a = 0.06, b = 0.4. Red line: f{t). Blue line: g{t). Green line: r{t) 
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Figure 7: a = 0.01, b = 0.4. Red line: f{t). Blue line: g{t). Green line: r{t) 

These results may appear surprising, but can be explained by looking at the structure 
of the states P/j^it) and their entanglement characterization. Indeed, one knows that the 
so-called entanglement of formation j^] is a monotonically increasing function of the 
concurrence |j2l| and, moreover, that it cannot increase under local quantum operations. 
This seems to conflict with the fact that the local action of 7t ® id may increase the 
concurrence and hence the entanglement of formation. However, there is no contradic- 
tion: in fact, 7t is not positive, whereas a quantum operation is by definition completely 



positive. It is of some interest to inspect in more detail how non-positivity may gener- 
ate entanglement by acting locally. Let us then consider the explicit expression of the 
entanglement of formation, 



where S{pl) is the von Neumann entropy of the marginal states of party 1 resulting from 
partial trace over party 2, p] :— Tr2(pi), obtained from those that contribute to the 
convex decomposition p = YliPiPi^ Pi > ^, J2iPi ^ ^■ 

Essentially, E{G (g) id[p]) < E{p) under a quantum operation G because any optimal 
decomposition 'Y^^PiPi of p achieving E{p) provides a decomposition 'Y^^Pi'Gi (g) id[pj] of 
G ® id[p], whence the minimum can only decrease. In this section wc have proved that 
there are times < s < i such that Pfx{s) — ^s® and = ^t® id[-P^] are well- 

defined entangled isotropic states satisfying E{P^{t)) > E{P^{s)). Prom the argument 
of above, such inequality can hold only if 7^ <S) id acting on an optimal decomposition of 
P^{s) does not provide a decomposition of -P^(i), namely only if at least one of the states 
optimally decomposing P^(s) are not mapped into density matrices by jt ® id. 

To summarize, by resorting to the slippage mechanism, one may force a non-positive 
time-evolution 74(8) id to map a class of entangled states into states at all times, but it may 
happen that this local action increases their entanglement because it is the positivity of 
the spectrum of other entangled states, outside that class, which is spoiled in the course 
of time. 



Wc have considered a concrete model of reduced dynamics for a single qubit weakly 
interacting with a stochastic environment and derived a reduced dynamics depending on 
two parameters. The resulting semigroup consists of maps 7^ that range from non-positive 
to positive and completely positive. 

The physical meaning of complete positivity is related to the existence of entangled 
states and can thus be fully appreciated only when the open qubit is correlated to an 
ancilla. Far from being abstract and out of experimental control, an entangled qubit pair 
one of whose parties experiences a noisy channel is a typical theoretical and practical 
context of quantum information and communication theory. In these cases, complete 




(51) 



VI. CONCLUSION 



positivity cannot be dispensed with, otherwise physical inconsistencies immediately ap- 
pear, typically the presence of negative probabilities in the spectrum of evolving entangled 
density matrices which lose their meaning as physical states. 

In relation to entanglement, we have studied an approach, the so-called slippage of 
initial conditions, whereby complete positivity is deemed an abstract request devoid of 
physical content and as such refused as an unnecessary constraint on the open quantum 
dynamics. This means that one prefers to stick to a non-positive Markovian reduced 
dynamics and try to cure in some way or the other its pathological behavior. Indeed, 
already for a single qubit, one has to enforce a selection of the admissible initial states 
which is supposedly due to the action of the transient regime. 

By means of a simple (completely positive) slippage operator S^, < /i < 1, we have 
turned non-positive reduced dynamical maps 74 into completely positive ones, 7f o 
that eliminate the presence of negative eigenvalues from the spectrum of time-evolving 
single qubit states. 

However, we have showed that 74 ® id acting on the isotropic states = ® id[P], 
may increase their entanglement by acting locally. Such an unphysical possibility is due 
to the non-positivity of 7t and can be eliminated by resorting to further slippage operators 
that may result in the whole elimination of entangled isotropic states. 

Despite the simplicity of the model, we believe it contains the salient features of a more 
general structure: sticking to non-completely positive reduced dynamics, though cured 
by some suitable slippage mechanism, would conflict with the presence of entanglement. 
In some particular cases discussed in this paper, the conflict can only be avoided by the 
drastic elimination, via a suitable slippage, of whole classes of entangled states. 

To conclude, slipped no n- (completely) positive semigroups present two different con- 
tradictory aspects in relation to entanglement: on one hand, the slippage operated by 
the bath during the transient phase seems to hamper the possibility of creating entan- 
glement by embedding two dynamically uncorrelated parties within a same environment, 
as discussed in On the other hand, unless all entangled states are eliminated by 
the slippage beforehand, the subsequent Markovian time-evolution may look like being 
able to create entanglement by acting locally. However, this entanglement creation is 
a spurious artifact: it is due to some entangled state, not destroyed by the slippage, 
which develops negative probabilities while evolving in time. As a consequence, curing 
non-positive Markovian evolutions through the slippage operation appears to be rather 



unsatisfactory and physically unviable. 
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